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We study the electronic structure of cubic phase uranium dioxide at different volumes using a 
combination of density functional theory and dynamical mean-field theory. The ab initio calculations 
predict an orbital-selective Mott insulator-metal transition at a moderate pressure of ~ 45 GPa. At 
this pressure the j = 5/2 states become metallic, while the j = 7/2 states remain insulating up to 
about 60 GPa. In the metallic state, we observe a rapid decrease of the 5/ occupation and total 
angular momentum with pressure. Simultaneously, the generalized Zhang-Rice state, which is of 
predominantly j — 5/2 character, quickly disappears after the transition into the metallic phase. 

PACS numbers: 71.15.-m, 71.20.-b, 71.27.+a, 71.30.+h 


Over the past decades, actinide materials including 
pure elements, hydrides, oxides, carbides, and nitrides 
have been extensively studied with numerous experimen¬ 
tal and theoretical tools, due to their fundamental impor¬ 
tance in the nuclear energy industry and military technol¬ 
ogy [1, 2]. Even though great progress has been made, 
many problems and puzzles still remain. Of particular 
interest is the electronic structure of these actinide ma¬ 
terials under extreme conditions (for example high pres¬ 
sure and high temperature in a reactor environment or a 
reactor accident) [3, 4]. 

Experimentally, it has been observed that under pres¬ 
sure many actinide materials undergo a series of struc¬ 
tural phase transitions, such as the cubic to orthorhombic 
phase transitions occurring in some actinide dioxides [5] , 
or the three successive phase transitions in Am [6]. To 
explain these complex phase transitions requires an ac¬ 
curate description of the properties of the 5/ electrons 
over a wide range of pressures. It is well known that the 
5/ electrons, which play a pivotal role in determining 
the key physical and chemical properties of the actinide 
materials, react sensitively to changes in the surrounding 
environment [1, 2]. So, it is natural to expect that the 
5/ electronic structure of actinide materials will be mod¬ 
ified and some exotic correlation phenomena may emerge 
when an external pressure is applied. Up to now, only a 
few experiments and calculations have been conducted to 
explore the high pressure properties of actinide materi¬ 
als, and most of these efforts were devoted to study their 
structural phase transitions and instabilities [4-7]. As 
a consequence, the high pressure electronic structures of 
actinide materials are poorly understood. Here, we em¬ 
ploy a state-of-the-art first-principles approach to shed 
new light onto the pressure-driven electronic transitions 
in uranium dioxide. 

Stoichiometric UO 2 is in a cubic fluorite structure at 
ambient pressure. A transition to the orthorhombic co- 


tunnite structure in the range of 42 ^ 69 GPa [5] was 
determined by high precision X-ray diffraction experi¬ 
ments. The ground state of UO 2 is an antiferromagnetic 
Mott insulator with a sizable band gap of ^ 2.1 eV [8]. 
To study the electronic structure of UO 2 , many tradi¬ 
tional first-principles approaches have been employed, 
such as the local density approximation (LDA) [or gener¬ 
alized gradient approximation (GGA)] plus Hubbard U 
approach [9-13], the hybrid functional method [14], and 
the self-interaction corrected local spin-density approx¬ 
imation (SIC-LSDA) [15]. However, none of the above 
methods can provide a satisfactory description of UO 2 
over a wide range of conditions. Eor instance, though the 
band gap at ambient pressure is reproduced correctly, the 
phase transition pressures predicted by the LDA + U and 
GGA + U methods are 7.8 [9] and 20 GPa [13], respec¬ 
tively, and thus much lower than the experimental value. 
This discrepancy is likely due to the incorrect treatment 
of the multiple metastable states of UO 2 [11]. In addi¬ 
tion, the changes of the electronic structure, especially 
the evolution of the Mott gap of UO 2 under pressure, is 
still an open and debated question [3]. 

In this letter, we use the density functional the¬ 
ory plus dynamical mean-field theory (DET + DMET) 
method [16, 17], which is currently the most powerful 
established tool to study strongly correlated systems, in¬ 
cluding the actinide materials [18-21], to calculate the 
electronic structure of cubic phase UO 2 under different 
pressures [22]. We predict that an orbital-selective Mott 
insulator-metal transition will happen around 45 GPa. 
We also find that the pressure effect will lead to not only 
a redistribution of the 5/ electrons and an associated 
collapse of total angular momentum, but also a disap¬ 
pearance of the generalized Zhang-Rice state (ZRS). 

Physical properties at ambient pressure. We first cal¬ 
culate the bulk properties of cubic phase UO 2 at am¬ 
bient pressure and T « 116 K [22] to benchmark the 
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FIG. 1. (color online). Electronic structure of cubic UO 2 at ambient pressure, (a) Total and 5/ partial density of states obtained 
by the DFT + SOC method, (b) Total and 5/ partial density of states obtained by the DFT + DMFT method {T ^ 116 K). 
The XPS and BIS experimental data (denoted by red hlled circles) [8] are shown for comparison, (c) The momentum-resolved 
spectral function obtained by the DFT + DMFT method. Here “ZRS” indicates the generalized Zhang-Rice state. 


TABLE I. Calculated bulk properties and band gaps of cu¬ 
bic UO 2 . Here ao denotes the lattice constant, B the bulk 
modulus, and B the hrst pressure derivative of B. 



ao (A) 

B (GPa) 

b' (GPa) 

gap 

. (eV) 

The present work 

5.571 


183.8 


4.24 

2.1 


LDA (GGA) + U 

5.540 

[10] 

191.6 

[10] 

- 

2.2 

[10] 


5.449 

[12] 

222.4 

[12] 

- 

2.3 

[12] 

Hybrid functional 

5.463 

[14] 

218.0 

[14] 

- 

2.4 

[14] 

SIG-LSDA 

5.400 

[15] 

219.0 

[15] 

- 

0.0 

[15] 

Experiments 

5.473 

[5] 

207.2 

[5] 

4.54 [5] 

2.1 

[8] 


accuracy of the DFT + DMFT approach. We use the 
equation of states (EOS) proposed by Teter et al [23] 
to fit the E-V curve, and then extract useful parameters. 
The calculated bulk properties are summarized in Tab. I. 
Other theoretical and experimental values, if available, 
are collected and listed as well. The obtained equilib¬ 
rium lattice parameter ao is slightly overestimated, while 
the bulk modulus B is underestimated. The small devi¬ 
ations (1.7% for ao, 11.3% for 5, and 6.7% for B ) may 
be caused by the large Hubbard U parameter used in the 
present calculations [22]. However, utilizing the obtained 
EOS, we find that at the experimental lattice volume, the 
calculated bulk modulus is 207.7 GPa, which is almost 
identical with the experimental one [5]. 

The electronic structure of UO 2 at ambient pressure is 
shown in Fig. 1. We can see that the DFT + spin-orbit 
couping (SOC) method fails to reproduce the insulat¬ 
ing nature of UO 2 , and leads to a metallic state. Thus, 
considering the SOC effect alone is not sufficient for a 


proper description of the electronic structure of UO 2 . In 
Fig. 1(b) and (c), the total (partial) density of states and 
momentum-resolved spectral functions obtained by the 
DFT + DMFT approach are shown, together with the 
available XPS and BIS experimental data [ 8 ]. From the 
calculated results, we can conclude that: (i) the energy 
gap is about 2.1 eV, which is consistent with the exper¬ 
imental value (see Tab. I) [ 8 ]. (ii) Since the energy gap 
is mainly associated with 5/ ^ 5/ transitions, uranium 
dioxide is a Mott insulator, which agrees with the experi¬ 
ments [24, 25] and most of the first-principles results [14]. 
(iii) Between —2.0 and —1.0 eV, there exists an isolated 
peak which is composed predominantly of U-5/ and O- 
2p characters. The position and spectral weight of it 
agree quite well with the XPS experiments [ 8 , 24]. Re¬ 
cently, Yin et al. [18] suggested that this low-energy reso¬ 
nance can be viewed as a generalized ZRS [26] . Since the 
many-body ground state of the UO 2 5/^ electronic con¬ 
figuration is a F 5 triplet, there must be a local magnetic 
moment which will couple with the 0 - 2 p’s hole induced 
by the photoemission process. Thus, in the following we 
identify this peak as the ZRS. 

Mott insulator-metal transition. Next, we concentrate 
on the pressure-driven electronic structure transition in 
cubic phase UO 2 . We decrease the lattice constant ao to 
mimic the increase of external pressure, while the crystal 
structure (including the space group symmetry, which is 
Fm3m) and the atomic coordinates are kept fixed. The 
evolution of the density of states with volume, as pre¬ 
dicted by DFT + DMFT, is shown in Fig. 2(a). As 
the volume (pressure) is decreased (increased), dramatic 
changes are observed in the density of states. The band 
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FIG. 2. (color online). Evolution of the electronic structure of cubic UO 2 under pressure calculated by the DFT + DMFT 
method, (a) Total density of states, (b) The U-5/ j = 5/2 partial density of states, (c) The U-5/ j = 7/2 partial density of 
states, (d)-(h) The momentum-resolved spectral functions for V/Vq — 1.00, 0.90, 0.85, 0.83, and 0.78, respectively. Notice the 
sharp quasi-particle peak pinned at the Fermi level in panel (f). 


gap shrinks monotonously, and disappears suddenly at 
V/Vq = 0.85, which marks a Mott insulator-metal tran¬ 
sition. The critical transition pressure Pc is about 45 
GPa according to the experimental P-V curve [5]. When 
P > Pc ^ the quasi-particle peak grows quickly, which 
implies an enhancement of metallicity. The peak associ¬ 
ated with the ZRS is also very sensitive to the volume 
collapse. Even though no strong shift of the ZRS peak 
to lower energies is apparent, which roughly agrees with 
previous DFT + DMFT calculations [18], there is a sub¬ 
stantial transfer of spectral weight to the quasi-particle 
peak. At V/Vq = 0.80 (corresponding to P ^ 65 GPa) [5], 
the ZRS peak is almost smeared out. The multiple peaks 
from —12 eV to —3 eV and from 1 eV to 10 eV are pre¬ 
dominantly of 0-2p character [24]. They are strongly 
hybridized with the lower and upper Hubbard bands of 
the U-5/ orbitals which are approximately located at —4 
to —5 eV and 1 to 6 eV, respectively [18]. These lig¬ 
and bands are also shifted outward and broaden under 


pressure. We further calculate the momentum-dependent 
spectral functions at various volumes in order to gain a 
better understanding of the evolution of the Mott gap 
and ZRS. Selected results are visualized in Fig. 2(d)-(h). 
These plots reveal the same Mott transition scenario as 
already shown in Fig. 2(a). From (d) to (f), the Mott gap 
is reduced from a finite value (2.1 eV) to zero, and the 
ZRS feature is slightly shifted toward the Fermi level. In 
panel (f), where U/Vb = 0.85, the quasi-particle peak ap¬ 
pears, while the ZRS peak still persists. From (f) to (h), 
the ZRS peak is rapidly smeared out and at the same 
time, the weight of the quasi-particle peak is substan¬ 
tially enhanced. 

Due to the SOC effect in heavy elements, the U-5/ or¬ 
bitals are split into two components: j = 5/2 and j = 7/2 
states (here we ignore the crystal field splitting for the 
sake of simplicity) [1]. It is interesting to study the evolu¬ 
tion of the corresponding density of states under pressure. 
The calculated results are shown in Fig. 2(b) and (c), re- 
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FIG. 3. (color online), (a) The average 5/ electron occupancy Nf = (N) (left y-axis) and total angular momentum (J) (right 
y-axis) for various volumes, (b)-(c) Valence state histograms at V/Vq = 1.00 (insulating region) and 0.78 (metallic region), 
respectively. The pie diagrams show the distributions of atomic states with respect to N. The atomic state probabilities for 
V = 0 and V = 4 states are too small to be visible. 


spectively. If one compares the j = 5/2 partial density of 
states in Fig. 2(b) with the corresponding total density 
of states in Fig. 2(a), one can easily recognize that both 
the Mott gap and the ZRS peak are associated primarily 
with the j = 5/2 states. In Fig. 2(b), on can observe 
again how a quasi-particle peak grows at the Fermi level 
and the ZRS peak fades away when V/Vb ^ 0.85. Thus, 
we propose that the j = 5/2 states are responsible for 
the pressure-driven Mott insulator-metal transition and 
the associated disappearance of the ZRS peak in cubic 
UO 2 . As for the j = 7/2 states, most of the density of 
states is well above the Fermi level. Remarkably, even 
when V/Vb = 0.85 (the approximate transition point for 
the j = 5/2 states), the spectral weight for the j = 7/2 
states at the Fermi level remains zero. Only for V/Vb < 
0.83, the j = 7/2 states become metallic. In other words, 
the Mott insulator-metal transitions for the j = 5/2 and 
j = 7/2 states don’t occur simultaneously, and there is 
a sizable volume or pressure range (0.83 < V/Vq < 0.85, 
corresponding to 45 GPa < P < 60 GPa [5]) in which 
the j = 5/2 states are metallic while the j = 7/2 states 
remain insulating. According to our calculations, the 
insulator-metal transition in cubic phase UO 2 is there¬ 
fore an orbital-selective Mott transition [27, 28]. 

Finally, we would like to emphasize that the structural 
phase transition is completely ignored in our calculations. 
In fact, the structural phase transition from cubic to or¬ 
thorhombic begins at V/Vb = 0.87, and continues beyond 
V/Vb = 0.81. The transition zone extends at least up to 
69 GPa [5]. Since the Mott transition apparently oc¬ 
curs without any change in structure and the Pc for it 
is very close to the one for the structural phase transi¬ 
tion, we believe that the structural phase transition in 


UO 2 is driven by the Mott transition, which involves a 
localization-delocalization process of the U-5/ electrons. 
Previous theoretical research using the LDA + U method 
predicted that the orthorhombic structure is an insula¬ 
tor, and the metallization pressure for it is in the range 
of 226 ^ 294 GPa [9]. On the other hand, the SIC-LSDA 
approach predicts that the cubic UO 2 is on the verge 
of an insulator-metal transition at ambient pressure [3]. 
Both predictions are inconsistent with ours and it would 
thus be very interesting to test them experimentally. 

5/ orbital occupancies and atomic multiplets. The 
nominal 5/ occupation for UO 2 is 2. However, due to 
the strong hybridization between the U-5/ and 0-2p 
bands [24], the actual 5/ occupation is larger than 2. 
We analyzed the connection between the average 5/ or¬ 
bital occupancy Nf and the volume collapse V/Vq [see 
Fig. 3(a)]. Under compression, Nf first increases steadily, 
and then decreases sharply. This behavior can be easily 
understood as follows: In the insulating phase, the U-5/ 
electrons are localized and the moderate changes in Nf 
are due to the hybridization effect with the 0-2p bands. 
In the metallic phase, the U-5/ electrons become itiner¬ 
ant which allows them to hop from the j = 5/2 states 
to the j = 7/2 states and the 0-2p orbitals. As a re¬ 
sult, one observes an abrupt decline of Nf after the Mott 
insulator-metal transition. 

The valence state histogram, which measures the prob¬ 
ability to find a 5/ electron in a given atomic eigen¬ 
state, can provide additional information about the na¬ 
ture of the U-5/ electrons. Representative valence state 
histograms for U/Vb = 1.00 (insulating phase) and 0.78 
(metallic phase) are plotted in Fig. 3(b) and (c), respec¬ 
tively. In the insulating phase, the dominated states have 
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occupancy N = 2 and 3. The other states have only 
negligible weight. In the metallic phase, even though 
the N = 2 and 3 states remain dominant, the probabili¬ 
ties for the N = 1 states increase substantially. So, the 
change in the valence state histogram is consistent with 
the rapid decrease of Nf . We further calculated the effec¬ 
tive total angular momentum (J) using (J) = 
where pr is the probability for the atomic eigenstate T. 
In the insulating phase, (J) is hardly affected by the vol¬ 
ume compression, while in the metallic phase it decreases 
steeply [see Fig. 3(a)]. We note that the collapse of the 
total angular momentum under pressure is rather similar 
to the spin state crossover observed in some transition 
metal oxides [29], so it is likely the driving force behind 
the “melting” of the ZRS peak [22] . 

In summary, we predict that the Mott insulator-metal 
transition in cubic UO2 occurs at Pc ~ 45 GPa and is 
of the orbital-selective type [27, 28]. There are two suc¬ 
cessive transitions for the j = 5/2 and j = 7/2 states, 
respectively. The ZRS peak is only prominent in the insu¬ 
lating phase. In the metallic state, it quickly fades away 
with increasing pressure. At the same time, in the metal¬ 
lic phase, the pressure effect leads to a rapid decrease of 
the 5/ occupancies and a collapse of the total angular 
moment. To our knowledge, this is the first study of the 
high pressure electronic structure of actinide dioxide us¬ 
ing a modern DFT -f DMFT approach. Our findings 
imply that physical properties such as resistivity, optical 
conductivity, and the magnetic moment of UO2 should 
change radically across the (electronically driven) struc¬ 
tural phase transition. We further speculate that other 
actinide dioxides, such as Pa02, Np02, PUO2, Am02, 
and Cm02 etc., should exhibit similar pressure-driven 
phenomena [14]. Thus, more high pressure experiments 
and theoretical calculations are highly desired. 
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